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Topological insulator represents a new class of quantum phase defined 
by invariant symmetries and spin-orbit coupling that guarantees metal- 
lic Dirac excitations at its surface. The discoveries of these states have 
sparked the hope of realizing nontrivial excitations and novel effects like 
magnetoelectric effect and topological Majorana excitations. Here we 
develop a theoretical formalism to show that three dimensional topo- 
logical insulator can be designed artificially via staking bilayers of two- 
dimensional Fermi gases with opposite Rashba-type spin-orbit coupling on 
adjacent layers, and with inter-layer quantum tunneling. We demonstrate 
that in the stack of bilayers grown along (OOl)-direction a nontrivial topo- 
logical phase transition occurs above a critical number of Rashba-bilayers. 
In the topological phase we find the formation of a single spin-polarized 
Dirac cone at the F- momentum point. This approach offers an accessible 
way to design artificial topological insulators in a set up that takes full 
advantage of atomic layer deposition approach. This approach also allows 
us to bypass limitations imposed by bulk crystal geometry and is tunable. 
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The unusual properties of topological insulator (TI) rely on searching suitable materials 
that inherit invariant symmetries and spin-orbit coupling in the bulk ground state, and 
that allow the formation of metallic Dirac fermions at the boundary. [l]-[6] From theoretical 
standpoint, the bulk symmetries include combination of low- lying odd parity orbit als, and 
time-reversal symmetry with a crystal geometry that can conspire an inverted band structure 
via strong spin-orbit coupling,^ [8j or sometimes crystal mirror symmetry. [3j The main 
focus of the search of TIs is thus limited on synthesis of bulk materials with these inherent 
characteristics. Compounds of choice are the ones (Bi-based compounds P, [10] and their 
potential functional variants P [TTHT3] ) with large bulk spin orbit insulating gap. The efforts 
have then been extended to manipulating the 'non-trivial' topological phase by driving a 
trivial topological system through topological phase transition via chemical doping. [T4l [T5] 
or by tuning lattice constant, [13j or by introducing broken symmetry quantum phases. [161 
\T7\ These materials offer tremendous opportunity for application as well as realization of 
numerous nontrivial properties such as anti-localization, unusual magneto-electric effect, 
fractional quantum statistics, and controlled electronic mass by applied magnetic field, or 
proximity to quantum orders. p!l [2l [T8142T] Yet, as exciting as these possibilities are, exposing 
the Dirac cone to the bulk insulating gap or even sometimes obtaining true bulk insulator 
have remained uncharted, and thus pose an active topic of intense research. 

We propose here an alternative approach to design TIs by combining a set of layers of two 
dimensional Fermi gases (2DFGs) with Rashba-type spin-orbit coupling. The surface spin- 
orbit locked metallic state, in the presence of interlayer quantum tunneling, translates into a 
bulk insulator with Z2-invariant topological properties, with Dirac excitations on the surface. 
The idea is to grow two counter-helical Rashba-planes ('Rashba-bilayer') — which hitherto 
imposes time-reversal invariance — along the (OOl)-axis with an interlayer distance that 
enables single-electron hopping between them. With an effective Hamiltonian, we observe 
that above a critical number of Rashba-bilayers, about 5-6 layers for a realistic parameter 
choice, the non-trivial TI phase emerges. The resulting gapless single-Dirac cone has a 
linear slope determined by the bulk Rashba-coupling strength, and thus tunable internally 
and externally. With the rapid expansion of surface growth techniques like molecular beam 
epitaxy, one would be able to explore a large set of compounds that might not be accessible 
in the bulk phase, and yet possess TI behaviors. We find that this 'layer by layer' approach 
of Rashba-bilayers has all the same known set of properties as the bulk TI; for example we 
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present a direct calculation of the Chern index in the bulk to support the topological nature 
of the resultant state. With the rapid expansion of surface growth techniques like molecular 
beam epitaxy, one would be able to explore a large set of compounds that might not be 
accessible in the bulk phase, and yet possess TI behaviors. 

In what follows, we lay a general framework for generating and manipulating 'home- 
made' three dimensional TI. The presence of Dirac cone at the Fermi level without the 
intervention of bulk state will also make it a attractive candidate to realize the fractional 
statistics and non-Abelions.fT^ In principle one can envision to use the approach we pro- 
pose and take it few steps further, e.g. it allows one to introduce layers of magnetic. ^20l [2T] 
superconducting!^ [191 12] other many body orders||23j within the topological matrix. 
Taken together, the search for TI materials which obey several symmetry properties and 
inherits spin-orbit coupling can thus be replaced with 'homemade' systems by generating 
spin-orbit coupling via external electron field, and imposing symmetry properties via manip- 
ulating heterostructure geometry which will be free from any particular crystal geometry. [2^ 

Rashba-bilayer heterostructure 

We start with depicting our basic idea in Fig. [T| and constructing the low-energy eflFective 
Hamiltonian. The cornerstone of this framework relies on gluing two counter-propagating 
Rashba-type spin-orbit coupled 2DFGs, denoted by ±a{k)^ close to each other such that 
finite quantum tunneling [i^(k)] is activated between them, as illustrated in Fig. [T^. Such 
Rashba-bilayer can easily be manufactured by creating potential gradient between two 
2DFGs with the help of gating, or by inserting polarized medium between them, among 
others. As opposed to a metallic single-Rashba 2DFG, the Rashba-bilayer governs an insu- 
lating gap, determined by the value of D(k)^ with a minimum gap at the F-momentum point. 
Then Rashba-bilayers can be grown along (001) direction with an inter-bilayer electron hop- 
ping, t^, which is required to be diflFerent than D{0) to eliminate the degeneracy in the band 
structure. Above a critical number of the Rashba-bilayers in such heterostructure setup, a 
bulk topological phase transition commences, which turns the quantum well band curvature 
into an inverted band dispersion or 'dent' structure in the valence Fermi sea, as illustrated 
in Fig. [T]l. Such inverted band dispersion is also observed in first-principles band structure 
calculations, [25] and angle-resolved photoemission spectroscopy (ARPES) measurements [H] 
in 3D bulk materials. Both analytically and numerically, we investigate below a widespread 
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realistic parameter space in which a non-trivial topological phase transition is governed with 
single Dirac cone, carrying all salient topological properties that are derived and realized 
earlier in bulk three-dimensional systems. [U [2j A detailed progression of the resulting band 
structure is given in the supplementary information (SI). 

Effective model of Rashba-bilayer heterostructure 

Based on the above-proposed setup, we now derive an effective low-energy theory for the 
origin of TI. For each 2DFG planes, an electron with momentum k experiences an effective 
'anisotropic' magnetic field, induced by an electric field E^^ which couples to its spin cr to 
give rise to Rashba-type spin-splitting electronic bands = /2rrf ±aR{k x cr), where aR 
is Rashba-coupling strength, controlled by external or internal electric field, and m* is the 
effective mass of electrons, k is defined in the 2D plane, and correspondingly cr stands for 
Pauli matrices in the spin-subspace. We take two such counter-propagating helical 2DFGs, 
Kr and /i^, close to each other, either on the same plane or out of plane, such that finite 
quantum tunneling D{k) between them becomes active. Then this setup is grown along 
z— axis with a non-spin- flip electron hopping between them, t^. Therefore, the general form 
of the effective Hamiltonian for number of semi-infinite layers within open-boundary 
condition can be expressed as 
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Each term in the above Hamiltonian is a 2 x 2 matrix. To keep the formalism general 
and readily tunable, we take 'anisotropic' hoppings between Hr and Hr as D{k) = {Dq + 
Mk^)l2x2^ which can also be thought of as Dirac mass and Newtonian mass, respectively, 
due to their fundamental impacts on the obtained band structure. The tunneling between 
two adjacent bilayers is T = tzhx2' is a 2 x 2 component zero vector and I2x2 is the 
identity matrix. Time-reversal invariance of the above Hamiltonian consequently emerges 
due to the fact that hR^—k) = hR^k)^ which is an essential criterion for the formation of 
helical edge state, endowing the system to Z2 classification. [21 [4]-[6l El |26] 

The values of parameters m* and aR depend on the type of Rashba-setup under consid- 
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eration, whereas those for the expansion parameters M and tz are controUable mainly by 
the details of the heterostructure. An elaborated band structure progression as a function of 
these parameters is given in SI. Here we discuss the emergence of topological bulk insulator 
with helical surface state for a representative value of m* = 0.1 eV~^A~^, and aR = 1.5 eVA. 
Such value and even higher value of Rashba parameter is readily achievable in Bi-based sur- 
face and quantum- well states, as well as in BiTel bulk systems. [3j The Dirac mass has the 
same meaning as used for 2D quantum spin Hall (QSH) insulator. [28] and 3D TIs,[8] in that 
it is responsible for opening a band gap between the two adjacent time-reversal terms 
and at the F-point inside the bulk setup. When several slabs of two 2DFGs are glued 
together with t^, the insulating gap persists in every layers except on the boundary. Two 
helical edge states develop at each boundary which meet at a gapless Dirac cone above a 
critical value of the number of layers, Nc. 

The value of Nc depends only on the Dirac mass term Dq and on interlayer tunneling 
but not on other terms, and thus controllable solely by the 'layer- by-layer' deposition 
procedure. We illustrate this case for a representative value of Dq = —25 meV, tz = —0.2 eV 
in Fig. [2j Above three or more layers of Rashba-bilayers, we observe, both numerically and 
analytically (see SI), that an inverted bulk insulating gap develops in the interior Rashba- 
bilayers, which reflects in the transformation of the quantum- well-like states of the inherent 
slabs into inverted dispersions resembling a 'dent' shape in the flUed Fermi sea. We recall 
that this 'dent '-like inverted valance band is a critical signature of non-trivial topological 
phase, and has been consistently evaluated via flrst-principles calculations, [25] and also been 
observed in 3D TI materials. [T4] 

To further quantify the 'strong' topological phase transition as a function of number of 
layers, we also compute the topological invariants u (or equivalent ly an axion angle param- 
eter 9) of the Hamiltonian in Eq. ([T]). Since the Hamiltonian is invariant under inversion 
symmetry, we can derive the essential parity operator for each pair of Kramers degener- 
ate valence bands from the constraint PH{k)P~^ = H{—k)^ with P = Gz ® ^2x2, for each 
Rashba-bilayers. Thus we evaluate the topological quantum index as[29j 

{-^y = X{X{{Mm), (2) 

i 1=1 

where i stands for the the time-reversal high-symmetry points on the Brillouin zone, and 
is the eigenvector for each layer /. For the parameter set, we indeed flnd that a parity 
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inversion occurs at A^c, which, according to Z2 topological criterion, [21 HI |8l [26] marks the 
emergence of topological insulating phase at this critical value of A^. Above this critical 
thickness, the parity value stabilies to the non-trivial TI.f30] 

A consequence of the Z2 TI is the presence of spin-polarized Dirac cone at the surface. 
Despite the emergence of 'non-trivial' topological phase above three bi-layers, a surface gap 
persists due to finite quantum tunneling between the two edges, and the massless Dirac 
point disappears above five layers for this parameter choice. Such formation of massive 
Dirac quasiparticles below a critical value of quintuple layers is observed in Bi2Se3 thin 
films,[22] and also in doped bulk TlBi(Si_^Se^)2|ni CS]. 

We tabulate the values of bulk and surface gaps as a function of various tuning parameters 
in Fig.[3| As expected the electron's mass m* and Rashba-coupling aji do not have significant 
effect on the gap values, and thus gives us an alternative approach for generating bulk TI 
beyond atomistic spin-orbit coupling and heavy electron mass. The Dirac mass D and 
interlayer tunneling provide the backbone for engineering the surface and bulk gaps, which 
are readily tunable via heterostructure details. In S'/, we show the corresponding energy 
dispersions. 

Surface Dirac cone properties 

We can formulate the Hamiltonian for the the edge state by using the theory of invariants. 
The time-reversal symmetry imposes that a spin-up state at momentum h must be entangled 
to a state with spin-down at —k via spin-orbit coupling. As a consequence of this, according 
to Fermion doubling theorem, gapless Dirac point is guaranteed at F-point. In what follows, 
two counter-propagating surface states form with one spin state from the upper Rashba-term 
(/i^), and an opposite spin state from the lower Rashba term (/i^), which are related by 
time-reversal operation. By projecting the full Hamiltonian in Eq. 1 onto a single slab, we 
can write down the surface Hamiltonian in this basis to the leading order in k as 

a surf = 0LR{(y''K - cr^K)- (3) 

Note that the similar surface dispersion is also obtained on 3D TI in Refs. [2, 8j, except 
that here its slope is solely determined by the Rashba-coupling strength. Therefore, the 
velocity of the Dirac fermions turns out to he v = au/U^ and thus externally tunable. For 
the value of Rashba-coupling constant as large as 3.8 eVA achieved to date in bulk BiTeI,[3] 
we get '^2::^5.8xl0^m/s, which is even larger than the value achieved to date in a 3D 
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TIs.[9l[3l] The deviation from the hnear slope of the surface state within the bulk insulating 
gap depends on the value of mass term D{0). For the parameter set used above, we obtain 
linear dispersion expanding about 0.3 eV on both sides of the Dirac point, which is also 
tunable (see SI). 

The spin-polarized surface states can be imaged directly by ARPES and tunneling spec- 
troscopies, as well as via transport measurements. An interesting transport parameter of 
Dirac fermions is the quantum Hall effect which can be achieved by inserting a magnetic 
layer in the heterostructure next to the boundary, so that it breaks time-reversal symmetry 
on the surface via proximity effect, but not in the bulk. [21 [20] We can express the surface 
and magnetic layer Hamiltonian in a combined 'd '-vector form as d • cr, where the time- 
reversal invariant components are dx = —oiRkx and dy = aRky^ and time-reversal breaking 
dz gives the Zeeman energy splitting due to induced magnetization. Without the pres- 
ence of a dz term, the winding number of the upper and lower chiral states, defined as 
C = ^ J dkx J dkyd • {dxd x dyd) (the notation 'hat' represents corresponding unit vec- 
tor) where the integration is performed in closed loop around the Dirac cone, cancel each 
others. However, when a broken time-reversal symmetry is imposed, C becomes equal to 
zba^, where is the component of the spin along the magnetic field orientation. In this 
case, QSH conductance becomes fractional in units of e'^/h as an = =F(|<^z|/2)(e^//i). For 
fully polarized spin-configuration along the magnetic field, the above equation generates 
half-integer anomalous quantum-Hall effect (QHE). The fractional or half-integer QHE is a 
trademark feature of Dirac systems, [U [21 [3ll [32] and can be used as a test of our proposal 
from transport perspective. An interesting consequence of the half-integer landau level is 
that as TR breaking mass term, dz^ approaches zero, the counter-propagating Landau levels 
move to zero energy, and the system can act as protected quantum Hall 'insulator' despite 
the presence of energy states at the Fermi energy [33j. This unique scenario has been dis- 
cussed theoretically for TIs, but never been realized experimentally due to the absence of 
isolated Dirac point at Ep^[31\ however, it can be realized in the present heterostructure 
setup. 

Outlook 

Our proposal clearly overcomes the limitation of searching suitable combinations of crystal 
geometry, and inherent wave function symmetries within a bulk system for the material 
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realization and device implementation of the topological properties. The present formalism 
is free from any particular crystal geometry, and thus provides a widespread playground for 
engineering 'homemade' TI and surface state properties. Another advantage of the present 
case is that here the edge state is exposed to the bulk insulating gap, instead of being buried 
inside the bulk Fermi sea, which has greatly limited the usage of existing 3D TIs. The het- 
erostructure, with depositing capability of one atomic layer at a time, can also accommodate 
magnetic and superconducting layers, beyond the conventional doping or proximity effects, 
and thus will promote a unambiguous detection of non-Abelian particles such as 'axions', 
Majorana fermions p!8l [201 [21 J without the distraction of vortex, impurity, or dopant states. 
Given that electron interaction is stronger in 2DFG,[23] novel broken symmetry phases are 
easy to yield inherently in its heterostructure than in a weakly correlated 3D TI. 

We are indebted to S. Basak and A. M. Black-Schaffer for numerous discussions. This 
work was supported, in part, by UCOP and by Los Alamos National Laboratory, of the US 
Department of Energy under Contract DE-AC52-06NA25396, VR and ERG, and benefited 
from the allocation of supercomputer time at NERSG. 
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FIG. 1. Rashba-bilayer of 2DFGs and its heterostructure setup, a, A bilayer combination 
of opposite Rashba-type spin-orbit coupling 2DFGs, denoted by — representing /i^ and 

/i^, respectively in the Hamiltonian in Eq. 1. D(k) gives the interlayer electron tunneling between 
them, b, As grown Rashba-bilayers with finite interlayer coupling, t^. c, d, Illustration of band 
dispersions for a bilayer 2DFG, and its heterostructure version, respectively. The emergence of 
inverted band curvature in the valence bulk band marks the topological phase transition (see text). 
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FIG. 2. Band progression and formation of gapless surface state, a-d, Evolution of band 
dispersion for single layer (IL) to six layers (6L) of Rashba-bilayer setup, e-f, Corresponding 3D 
view of the band dispersions in the kx^ fc^-plane for the bulk (top panel) and surface (bottom panel) 
bands, respectively. The parameter set for this calculation is given in text. In going from two to 
three bilayers, the bulk valence band topology reveals the emergence of an inverted shape [see also 
c and corresponding f], which indicates the topological phase transition from trivial to non-trivial 
phase. However, it takes about six bilayers to turn off the inter-edge tunneling to commence gapless 
Dirac cone at the surface. Arrow dictate the spin orientation. The definite spin-chirality for the 
gapless Dirac cone is illustrated by counter-rotating arrows in f. 
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FIG. 3. Bulk and surface gaps in a realistic parameter space, a, The critical number of 
layers, Nc, at which an inverted bulk gap opens. The bulk band gap jumps from D(0) at A^=l 
(equivalent to ^^=0) to a large positive value, and then consistently becomes negative for > 3, 
constrained by the Hamiltonian (see SI). For each parameter set, the corresponding other parame- 
ters are set to be constants to the values mentioned in the main text. The varying parameters are: 
m*=l (green), 2.5 (blue), 5 (red) in eV~-'^A~^. aR^l (orange), 1.5 (green), 2 (blue), and 3 (red) in 
eV A which are attainable in existing system. [3J Dq=-25 (orange), -50 (green), -75 (blue), and -100 
(red) in meV, and t^=-10 (green), -20 (blue) and -50 (red) meV. Since Mk^ = at the F-point, 
the parameter M does not have any effect on the direct gap structure, b. Corresponding surface 
gaps at the F-point. Clearly, the gapless surface Dirac cone formation depends on two parameters, 
the interlayer hybridization t^, and strongly on the Dirac mass D. The parameter values here are 
same as in a. 
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BAND DISPERSIONS FOR A RASHBA-BILAYERS 



In this supplementary information, we provide further insights into the emergence of TI 
from analytical description. The Rashba-type spin-orbit coupling splits a non-interacting 
band into two in-plane spin-polarized bands. To keep the formalism general, we take a 
two-dimension Fermi gas (2DFG), with quadratic band, /c^/2m*, where m* is the effective 
mass of electrons. Under Rashba-coupling,[lJ the split bands are obtained by solving h% = 
/c^/2m*/2x2 ± OLR{ky(Jx + kx(Ty)^ where aR is Rashba-coupling. For simplicity, we define the 
second term as a{k). Next we create a bi-layer 2DFG with two counter-propagating helical 
2DEGs, and in which finite quantum tunneling is turned on, as denoted by D{k). 
The Hamiltonian for bi-layer 2DFG can be written explicitly as 



H 



k'^/rrf a{k) 

a\k) k^/rrf 
D{k) 

D{k) 



D{k) 

D{k) 

k'^/rrf —a{k) 

-a\k) k'^/rrf 



\ 



(4) 



The resulting eigenvalues are two- fold degenerate: E^{k) = k^ /rrf ± ^Ja^{k) + D'^{k). 
This specific bi-layer setup helps generating two upward and downward dispersing bands, 
with a direct gap between them defined by D{k) at F-point (since a(0) = 0), see Fig. Sib. 
i^(k) shall be even under inversion, e.g. D{k) = + Mk^^ which destroys the spin- 
polarization, implying that the gap between the band bands is charge insulating. 
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Fig. SI. Band progression from single to bilayer to heterostructure of Rashba-type spin-orbit 

coupling 2DFGs. 



BAND DISPERSIONS FOR TWO RASHBA-BILAYERS 



As depicted in Fig. 1 of the main text, we now take two bilayer 2DEG, grown along (001) 
direction such that quantum tunneling, say t^, between the adjacent layers becomes active. 
The explicit form of the Hamiltonian for this case is 
/ 
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The eigenvalues of Eq. [5] is given by 



Et{k) 



1 



m* 2 



2tl + 4 [a^{k) + D\k)) ± 2t,^tl + AD''{k) 



1/2 



(5) 



(6) 



The energy dispersions imply that whereas a large bulk gap opens between and Ez , the 
gap between the lowest energy states E^ and E^^ decreases from their single bilayer value of 
D^. Although, mathematically, there exists a solution for the values of and tz at which 
surface gap can be closed, but there is no clear evidence of the bulk band party inversion 
for this Hamiltonian, and thus the resulting solution remains trivial topological phase. The 
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absence of band inversion is also evident in the parity calculation and the parabolic shape 
of band dispersion as presented in the main text. 

BAND DISPERSIONS FOR THREE RASHBA-BILAYERS 

With three or more bilayers, coupled by finite quantum tunneling, the bulk systems begin 
to form in which the interior bilayer(s) acts as bulk where the edge layers give rise to surface 
states. We show here for three bilayers, a bulk band inversion occurs and thus promotes 
non-trivial topological insulating phase. However, to get rid of the electron hopping from 
the two edges, or in other words, to close the surface gap, one requires a higher number of 
bilayers. 

The analytical solution gradually becomes complicated and long with further increase of 
bilayers. However we can gain some insights into how an inversion occurs in three bilayers 
of 2DFGs by the symmetry of the Hamiltonian, irrespective of parameter choices. The 
Hamiltonian can be easily generalized from Eq. [5} and we present the eigenvalue for some 
of the relevant bulk bands as, 

E^{k) = ^± ^ [3 {3a\k) + 2tl + 3D{k)) X\k) + {tl + QD{k)) X{k) 



+Q{?>Y{k)-2)tl + Mt\D'{k)\"\ (7) 



where 



r 1 1/3 

X{k) = t, [-8tl + m,D{k) + l2Y{k)\ , 
and Y{k) = D{k) [-96D^(fe) - 39tlD\k) - Ut^] . (8) 

The important message of Eq.[7|is that there is a prefactor 1/X ^ which enables an 

inverted band gap opening for a negative value of t^. The resulting parity inversion in the 
valence band endows a non-trivial topological phase transition at a large parameter space 
(see Fig. 3 of main text). A visual proof of the inverted band gap can be obtained from the 
associated inverted band curvature of 'dent' shape in the valence bulk band, as illustrated 
in Fig. Sic. According to Fu-Kane criterion, [2j for such case gapless Dirac states appear at 
the surface. However, for most values of the parameters shown in supplementary Fig. S2, 
finite gap at the surface still persists due to finite size-eflFect, and the surface gap closing 
requires characteristic number of layers of Rashba-bilayers. 



17 



DEPENDENCE OF BULK AND SURFACE BAND STRUCTURE ON VARIOUS 

PARAMETERS 
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Fig. SI. Variations of dispersion relation for a wide range of parameter space, a, The 

band dispersion for a bilayer 2DFG to six layers setup as shown in Fig. SI of the main text. In 

the rest of the figures, only a single parameter is tuned (in each row), while keeping others 
parameters same to that in a. In each figures, the horizontal axis spans within kx = -0.2 to 0.2 
with ky — {) and the vertical axis runs from E —-OA to 0.4 eV with Fermi level {Ep) lying at 
the center (black dashed line), except in f2, where the energy scale runs from -0.6 to 0.6 eV. The 

value of Nc for the formation of gapless surface state depends only on tz and D^^ whereas the 
slope of the surface state is controlled by the Rashba-coupling an (compare b, cl, and c2), and 
the span of the linearly dispersing surface state depends on the bulk insulating gap, as well as on 
the Newtonian mass M, as expected form Eq. 2 of the surface Hamiltonian. 



For the benefit of understanding the role of each parameter on the bulk and surface band 
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structures, and also for the purpose of material growth, we present results for systematic 
tuning of various parameters as a function of number of layers in Fig. S2. The same results 
present in Fig. 2 of the main text is reproduced here in Fig. S2a and the parameters are 
listed in the left side of this panel. For each row in the figures below, a single parameter is 
tuned while the others are kept to the value listed in Fig. S2a. The eflFective mass, m* of 
the electrons does have any significant eflFect neither on the bulk and surface band-structure, 
nor on the slope of the emergent surface state. Since this quadratic term goes to zero at 
the F-term, it begins to play a role at higher-energy, and with increasing m*, it only slightly 
reduces the indirect bulk gap. The Rashba-coupling strength a^^ as for m*, has no role 
to the surface and bulk gaps (very minor contribution to the indirect bulk gap). However, 
as deduced in the main text in Eq. 3, determines the velocity of the Dirac fermions on 
the surface. Since is easily tunable both externally and internally, as well as, a large 
value of it can nowadays be achieved in available materials, [3j obtaining spin-polarized Dirac 
excitations with larger velocity than the presently known bulk 3D TI is possible. 

The Dirac mass modifies the band gap as its definition implies, in that it controls 
the surface gap, however, it does not have significant eflFect on the bulk gap. The band 
curvature of the surface state also remains very much unaltered. The Newton mass M is 
associated with a quadratic momentum dependence, and thus does not participate in the 
gap, and it contributes at high-energy. Higher the value of M, better the band topology. 
The inter-bilayer tunneling tz contributes almost oppositely to that of Dirac mass D^^ in 
that it modulates the bulk gap, however, it does not any prominent effect on the surface 
gap. 
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